Decoherence in two-dimensional quantum walks using four- and two-state particles 
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We study the decoherence effects originating from state nipping and depolarization for two- 
dimensional discrete-time quantum walks using four- state and two-state particles. By quantifying 
the quantum correlations between the particle and position degree of freedom and between the two 
spatial (x — y) degrees of freedom using measurement induced disturbance (MID), we show that 
the two schemes using a two-state particle are more robust against decoherence than the Grover 
walk, which uses a four-state particle. We also show that the symmetries which hold for two-state 
quantum walks breakdown for the Grover walk, adding to the various other advantages of using 
two-state particles over four- state particles. 



I. INTRODUCTION 

Quantum walks are a close quantum analog of classi- 
cal random walks, in which the evolution of a particle is 
given by a series of superpositions in position space p]- 
[5|. Recently they have emerged as an efficient tool to 
carry out quantum algorithms [6j [7] and have been sug- 
gested as an explanation for wavelike energy transfer 
within photosynthetic systems [5J [9] . They have applica- 
tions in the coherent control of atoms and Bose-Einstein 
condensates in optical lattices [lOj [11] , the creation of 
topological phases [12J, and the generation of entangle- 
ment [13 j. Quantum walks therefore have the potential 
to serve as a framework to simulate, control and under- 
stand the dynamics of a variety of physical and biologi- 
cal systems. Experimental implementations of quantum 
walks in last few years have included NMR [T4UT6] , cold 
ions p~7j [18] , photons [TM24] , and ultracold atoms [25] , 
which has drawn further interest of the wider scientific 
community to their study. 

The two most commonly studied forms of quantum 
walks are the continuous-time [26] and the discrete-time 
evolutions [5 , 27-32 . In this work we will focus on the 
discrete- time quantum walk and just call it quantum walk 
for simplicity. If we consider a one-dimensional (ID) ex- 
ample of a two-state particle initially in the state 

|*in) = (cos(o/2)|0) + e 4 " sin(«V2)|l)) ® |Vo), (1) 



then the operators to implement the walk are denned 
on the coin (particle) Hilbert space 1~L C and the position 
Hilbert space H p [H = H c Hp] . A full step is given by 
first using the the unitary quantum coin 
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and then following it by a conditional shift operation 

s x = [|o> <oi ® \^ x -i)(ip x \ + ii)<ii ® i^+ixvfei] • 

X 

(3) 

The state after t steps of evolution is therefore given by 

|* t > = [4[B(fl)®l]] t |*in). (4) 

All experimental implementations of quantum walks re- 
ported by today have used effectively ID dynamics. A 
natural extension of ID quantum walks to higher dimen- 
sion is to enlarge the Hilbert space of the particle with 
one basis state for each possible direction of evolution at 
the vertices. Therefore, the evolution has to be defined 
using an enlarged coin operation followed by an enlarged 
conditioned shift operation. For a two-dimensional (2D) 
rectangular lattice the dimension of the Hilbert space of 
the particle will be four and a four dimensional coin op- 
eration has to be used. Two examples of this are given 
by using either the degree four discrete Fourier opera- 
tor (DFO) [Fourier walk] or the Grover diffusion opera- 
tor (GDO) [Grover walk] as coin operations [33-35 . An 
alternative extension to two and higher (d) dimensions 
is to use d coupled qubits as internal states to evolve 
the walk [36, 37 . Both these methods are experimentally 
demanding and beyond the capability of current exper- 
imental set ups. Surprisingly however, two alternative 
schemes to implement quantum walks on a 2D lattice 
were recently proposed which use only two-state parti- 
cles. In one of these a single two-state particle is evolved 
in one dimension followed by the evolution in other di- 
mension using a Hadamard coin operation [38j[39]. In the 
other, a two- state particle is evolved in one dimension fol- 
lowed by the evolution in the other using basis states of 
different Pauli operators as translational states [40j [41] . 

In this work we expand the understanding of 2D quan- 
tum walks by studying the effects decoherence has on 
the four-state Grover walk and the two two-state walks 
mentioned above. The environmental effects are mod- 
eled using a state-flip and a depolarizing channel and we 
quantify the quantum correlations using a measure based 
on the disturbance induced by local measurements [42] . 
While in the absence of noise the probability distribu- 
tions for all three schemes are identical, the quantum 
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correlations built up during the evolutions differ signifi- 
cantly. However, due to the difference in the size of the 
particles Hilbert space for the Grover walk and the two- 
state walks, quantum correlations generated between the 
particle and the position space cannot be compared. The 
quantum correlations between the two spatial dimensions 
(x — y), obtained after tracing out the particle state, on 
the other hand, can be compared and we will show that 
they are larger for the walks using the two-state parti- 
cles. When taking the environmental effects into account, 
we find that all three schemes lead to different proba- 
bility distribution and decoherence is strongest for the 
Grover walk, therefore making the two-state walks more 
robust for maintaining quantum correlations. Interest- 
ingly, we also find that certain symmetries which hold 
for the two-state quantum walk in the presence of noise 
do not hold for a Grover walk. Together with the spe- 
cific initial state and the coin operation required for the 
evolution of the Grover walk, this reduces the chances to 
identify an equivalence class of operations on a four-state 
particle to help experimentally implement the quantum 
walk in any physical system that allows to manipulate 
the four internal states of the coin. 

This article is organized as follows : In Section[ll| we 
define the three schemes for the 2D quantum walk used 



to study the decoherence and in Section [III| we define the 
measure we use to quantify the quantum correlations. 
In Section IV , the effect of decoherence in the presence 



of a state-flip noise channel and a depolarizing channel 
are presented and we compare the quantum correlations 
between the x and y directions for the three schemes. We 
finally show in Section |V|that the state-flip and phase- 
flip symmetries, which hold for the two-state quantum 
walk, breakdown for the four-state walk and conclude in 
SectionED 



II. TWO-DIMENSIONAL QUANTUM WALKS 

A. Grover walk 

For a Grover walk of degree four the coin operation is 
given bv [331 [34] 
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and the shift operator is 
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where \i/j x ,y) = \^x) ® \^y)- It is well known that the 
operation [Sr x ^ [G0l]] t results in maximal spread of the 
probability distribution only for the very specific initial 
state 

|*f n ) = i(|0)-|l)-|2) + |3))®|^ 0) o), (7) 

whereas the walk is localized at the origin for any other 
case [43] [44] . Choosing |\£f n ) and evolving it for t steps 
one finds 

t t 

i* 4 w) = E E[A^),tio>+%,y),*ii> 

x= — t y=—t 

+C ( ^ ) , £ |2)+P ( ^ ) , £ |3)]0|^ ) ) (8) 

where A( Xiy ) it , ^(x, y ),t, c (x, y ),t and V^ y ) it are given by 
the iterative relations 



A 



1 



(x,y),t 



\-A( x+ 



l,2/+l),t-l 



l,y+l),t-l 



B 



(x,y),t 



c 



(x,y),t 



+C(cc+1,2/+1),£-1 + (9a) 

2^(^+1^-1)^-1 _ &(x+l,y-l),t-l 

(x-l,y-l),t-l + ^(x-\-l,y-l),t-l] (^b) 

-[A( x -i iy +i) it -l + B( x -l,y+l),t-l 

] (9c) 



x-l,y-l),t-l 



B { 



(x-l,y-l),t-l 



+C(x 



(x-l,y-l),t-l 



(x-l,y-l),t-l\- 



(9d) 



This results in the probability distribution 



P ^= E E [\ A ^v)A 
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which is shown in Fig. [I] for t = 25. 
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B. Alternate walk 

Very recently a 2D quantum walk was suggested which 
used only a two-state particle which walks first only along 
the x-axis followed by a step along the ?/-axis [38]. This 
walk can results in the same probability distribution as 
the Grover walk and its evolution is given by 

l*i> = (S(o,y)[B(6) ® i])(5 (X)0) [B(6>) ® i])|* in ), (11) 
where 

x,y 

+ \l)(l\®\ip x+1 , y )(ip x j] (12a) 
S(o,v) = E [l°X°l® \^,y-i)(->Px,y\ 

x,y 

+ |l)(l|®|V'*,»+l><V'*,y|]. ( 12b ) 
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FIG. 1: (Color online) Probability distribution for the Grover 
after 25 steps. An identical distribution using the two-state 
particle with the ^^"^ can be obtained using the alternate 
walk with the coin operation B(tt/4) and the Pauli walk after 
25 steps. 



Using a coin operation with = 7r/4, the state of the 
walk after t steps can then be calculated as 

|tf(t)) = W(7T/4)*|* in ) 

t t 

= Yl + B (*,2/),*I 1 >] ® \*P(x,y)}, 

(13) 
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where W(n/4) = (S {0ty) [B(n/4)®l])(S {xfi) [B(n/4)®i]), 
and A( x ^ y \j and B( Xjy \t are given by the coupled iterative 
relations 
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The resulting probability distribution is then 

t t 
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which for the initial state |^ in ) = ^(|0) + z| 1>) <8> |^o,o) 
gives the same probability distribution as the four-state 
Grover walk (see Fig.[T]). 



C. Pauli walk 

A further scheme to implement a 2D quantum walk 
using only a two-state particle can be constructed using 
different Pauli basis states as translational states for the 
two axis. For convenience we can choose the eigenstates 

"1 0" 



of the Pauli operator 03 
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, |0) and |1) as basis 
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states for x— axis and eigenstates of &i 

^(|0) + |1)) and |-) = ^(|0) - |1)) as basis states 
for y— axis [40]. which also implies that |0) = ^(|+) + 
|— )) and |1) = -^(|+) — | — )). In this scheme a coin 
operation is not necessary and each step of the walk can 
be implemented by the operation 



S a3 — 5(^,0) 
followed by the operation 

S*. = > ' \\+)(+\®\^x,y-l)(lpx,y\ 
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The state after t steps of quantum walk is then given by 

\9 t ) = [S^sY^in) 
t t 

= E E [^,»),*i°)+ B («,»),ti i >]®i^,»). 
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where A( XyV ) yt and Br x ^ y ^ t are given by the coupled iter- 
ative relations 
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The probability distribution 
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(20) 



is again equivalent to the distribution obtained using the 
Grover walk and therefore also to the alternative walk 
for the initial state |^ in ) = ^(|0) + i| 1>) <g> |^o,o) (see 

Fig-Q. 

While the shift operator for the Grover walk is defined 
by a single operation, experimentally it has to be imple- 
mented as a two shift operations. For example, to shift 
the state |0) from (x,y) to (x — 1, y — 1), it has first to be 
shifted along one axis followed by the other, very simi- 
larly to the way it is done in the two-state quantum walk 
schemes. Therefore, a two-state quantum walk in 2D has 
many advantages over a four-state quantum walk. The 
two-state walk using different Pauli basis states for the 
different axes has the further advantage of not requir- 
ing a coin operation at all, making the experimental task 
even simpler in physical systems where access to differ- 
ent Pauli basis states as translational states is available. 
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One can, of course, also consider including a coin opera- 
tion in the Pauli walk, which would result in a different 
probability distributions [4T] . 

For general initial states, quantum walks in 2D with a 
coin operation G U (2) can result in a many non-localized 
probability distribution in position space. This is a fur- 
ther difference to the Grover walk which is very specific 
with respect to the initial state of the four-state particle 
and the coin operation. 



III. MEASUREMENT INDUCED 
DISTURBANCE 

Quantifying non-classical correlations inherent in a cer- 
tain state is currently one of the most actively studied 
topics in physics (see for example [45]). While many 
of the suggested methods involve optimization, making 
them computationally hard, Luo [42] recently proposed a 
computable measure that avoids this complication: if one 
considers a bipartite state p living in the Hilbert space 
T-La^Hb, one can define a reasonable measure of the to- 
tal correlations between the systems A and B using the 
mutual information 

Hp) = s( PA ) + s( PB ) - s( P ), (21) 

where S(-) denotes von Neumann entropy and pA and 
pB are the respective reduced density matrices. If pA = 
^2jP J A^A an d 9b — YljPB^-B' then the measurement 
induced by the spectral components of the reduced states 
is 

n(p) = 2ni®n|pni®n|. (22) 

Given that 7[II(p)] is a good measure of classical correla- 
tions in p, one may consider a measure for quantum cor- 
relations defined by the so-called Measurement Induced 
Disturbance (MID) |42] 

Q(p) = I(p) - I[U(p)}. (23) 

MID does not involve any optimization over local mea- 
surements and can be seen as a loose upper bound on 
quantum discord [46 . At the same time it is known to 
capture most of the detailed trends in the behaviour of 
quantum correlations during quantum walks [47] . There- 
fore, we will use MID (Q(p)) in the following to quantify 
quantum correlations for the different 2D quantum walk 
evolutions. 

Despite having the same probability distributions in 
the absence of noise, the MIDs for the four-state walk 
and the two-state walks differ. In Fig.|2|a) we show the 
MID between the particle and the position degree of free- 
dom, Q(pp P ) for all three walks and find that it is signif- 
icantly higher for the Grover walk. Where p pp is p4 S (t) 
for Grover walk, p2 S (t) for alternate walk and p2sa(t) f° r 
the Pauli walk. However, due to the difference in the 
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FIG. 2: (Color online) Quantum correlations Q(p PP ) and 
Q(pxy) for the different schemes in the absence of noise. 

size of the particles degree of freedom for the Grover and 
the two-state walks, a direct comparison of the quantum 
correlations Q(p pp ) does not make sense. Among the two- 
state schemes on the other hand, we see that the Pauli 
walk has a larger Q(p pp ) in comparison to the alternate 
walk. 

A fair comparison between all systems can be made 
by looking at the quantum correlations generated be- 
tween the two spatial dimensions x and y, Q(p xy ) (see 
Fig.J^b)). p xy is obtained by tracing out the particle 
degree of freedom from from complete density matrix 
[p4s(t), p2s(t) and p2 Sa (t)] comprising of the particle and 
the position space. We find that Q(p xy ) is identical for 
both two-state schemes and exceeding the Grover walk 
result. This behaviour is similar to the one described in 
Refs. [38, 39 , where the entanglement created during 
the Grover walk was compared with the alternate walk 
using the negativity of the partial transpose, in its gen- 
eralization for higher-dimensional systems [48j [49] . 

IV. DECOHERENCE 

The effects of noise on ID quantum walks has been 
widely studied [47[ I50H53] but the implications in 2D 
settings are less well known [37, 54, 55 . In particular no 
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(a) (b) 



FIG. 3: (Color online) Probability distribution of the Grover 
walk when subjected to a different state-flip noise level with 
k = 23 after 15 steps, (a) and (b) are for noise levels, p = 0.1 
and p = 0.9, respectively. 



study has been done on either of the two-state schemes 
presented in the previous section and we therefore now 
compare their decoherence properties to the Grover walk, 
using a state-flip and a depolarizing channel as noise 
models. We show that this leads to differing probability 
distributions and has an effect on the amounts of quan- 
tum correlations as well. 



A. State flip noise 
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1. Grover walk 

For a two-state particle, state-flip noise simply induces 
"0 I 



a bit flip [cri 



1 



but for the Grover walk, the state- 



flip noise on the four basis states can change one state to 
23 other possible permutations. Therefore, the density 
matrix after t steps in the presence of a state-flip noise 
channel can be written as 



P4s(t) 



Hi- P )[s 4 p 4s (t-i)s 4 % 



(24) 



where p is the noise level, 64 = S( Xjy \[G 1], and the fi 
are the state-flip operations. For a noisy chann el with all 
the 2 3 possible flips one has k = 23 and in Figs. |3(a)] and 
3(b) we show the probability distribution of the Grover 
for weak (p = 0.1) and strong (p = 0.9) noise levels after 
15 step. Compared to the distribution in the absence of 
noise (see Fig. [T]) a progressive reduction in the quantum 
spread is clearly visible. Note that for p = 1 the walk 
corresponds to a fully classical evolution. 

In Fig.|4^a) we show the quantum correlations between 
the particle state with the position space, Q(p pp ), as a 
function of number of steps t. With increasing noise level, 
a decrease in Q(p pp ) is seen, whereas for the quantum cor- 
relations between the x and y spatial dimensions, Q(p xy ), 
the same amount of noise mainly leads to a decrease in 
the positive slope (see Fig.Qb)). 



FIG. 4: (Color online) Quantum correlations created by the 
Grover walk in the presence of a noise channel including all 
possible state flips (k = 23). 



2. Two-state walks 

The evolution of each step of the two-state quantum 
walk comprises of a move along one axis followed a move 
along the other. Therefore, the walk can be subjected to 
a noise channel after evolution along each axis or after 
each full step of the walk. In the first case, the noise level 
p' = I is applied two times during each step in oder to be 
equivalent to the application of a noise of strength p in 
the second case. For the alternate walk the evolution of 
the density matrix with a bit-flip noise channel applied 
after evolution along each axis is then given by 



P2s(t) 



(l-f) [S x p 2s (t-l)S x 
(l-f) [SyP2s(t)Sl 
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(25b) 



where g\ 
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1, Sy = S(o ty) [B{6) (8)1], and S x 



S( XtO) [B(0)®i]. 

Similarly, the density matrix with a bit-flip noise ap- 
plied after the evolution along each axis for the Pauli 



(a) Alternate walk 



(b) Alternate walk 
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FIG. 5: (Color online) Probability distributions of the two- 
state walk with the bit-flip noise applied after evolution in 
each direction. The strength of the noise is p/2 — 0.05 in (a) 
and (c) and p/2 = 0.45 in (c) and (d) and the evolution was 
carried out for 25 step each time. 



walk is given by 
P 
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In Figs. |5 (a)| and [5(b)] the probability distributions for 
the alternate walk after 25 steps of noisy evolution with 



p = 0.1 and p = 0.9 are shown and Figs. 5(c) and 5(d) 
show the same for the Pauli walk. It can be seen that 
the bit-flip noise channel acts symmetrically on both axes 
for the alternate walk, but asymmetrically on the Pauli 
walk. This is due to the fact that the bit-flip noise applied 
along the axis in which the <j\ Pauli basis is used leaves 
the state unchanged. A completely classical evolution is 
recovered for p = 1 (p f = 0.5) and an evolution with 
p = 2 is equivalent to one with p = 0. 

Evolving the density matrix and calculating the MID 
for a noiseless evolution (p = 0), one can see from Fig. [6] 
that the initial difference in Q(p pp ) between the Pauli 
walk and the alternate walk decreases during the evolu- 
tion and eventually both values settle around 1.5. For a 
noisy evolution, however, the initial difference in Q{p pp ) 
does not decrease over time and we find a higher value 
for the Pauli walk compared to the alternate walk. Simi- 
larly, careful examination of Fig. [7] shows that the Q(p xy ) 
for the alternate walk and the Pauli walk are identical in 
the absence of noise, but differ for noisy evolution, with 
the alternate walk being affected stronger than the Pauli 
walk. 
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(b) Pauli Walk 

FIG. 6: (Color online) Particle-position quantum correlations 
created by the two-state walks for different bit-flip noise levels 
applied after evolution along each axis. 



The density matrix for the second case, that is, with 
a noisy channel applied only once after one full step of 
walk evolution, for both two-state walks is given by 



(l-p)Q(t-l) (27a) 

+" (27b) 



where 



Q(t - 1) = S y S x p 2s (t - l)StSl (28a) 
Q*(t - 1) = S ai S a3 p 2s a(t - 1)5+3^. (28b) 

In this case maximum decoherence and a completely clas- 
sical evolution is obtained for p = 0.5 and the evolution 
with p = 1 is equivalent to the one with p = 0. The ob- 
tained probability distributions are almost identical for 
both walks and differ only slightly from the ones obtained 
for the alternate walk with noise applied after evolution 
along each axis (see Fig. 5(a)| and |5(b)| ). The correlation 
functions Q(p pp ) and Q(p xy ) behave very similar for both 
walks (see Figs. [8] and [9| and we can conclude that the 
presence of bit-flip noise on both two-state walks when 
applied after a full step leads to equally strong decoher- 
ence. 
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FIG. 7: (Color online) Quantum correlations between the x 
and y spatial dimensions created by the two-state walks for 
different bit-flip noise levels applied after evolution along each 
axis 



FIG. 8: (Color online) Particle-position quantum correla- 
tions created by the two-state walks for different bit-flip noise 
levels applied after evolution of one complete step. 



B. Depolarizing channel 

To describe depolarizing noise we use the standard 
model in which the density matrix of our two-state sys- 
tem is replaced by a linear combination of a completely 
mixed and an unchanged state, 

P 

p = - (<7i/9<7i + <7 2 /9<72 + <7 3 /3<7 3 ) + (1 - p)p, (29) 

where <ti, &2 and 03 are the standard Pauli operators. To 
be able to compare the effects of the depolarizing channel 
on the Grover walk and the two-state walks we will apply 
the noise only once after each full step. 



1. Grover walk 

For the four-state particle the depolarizing noise chan- 
nel comprises of all possible state flips, phase flips and 
their combinations. State-flip noise alone leads to 23 
possible changes in the four-state system and adding the 
phase-flip noise and all combinations of these two is un- 
fortunately a task beyond current computational ability. 
Therefore let us first briefly investigate the possibility of 



approximating the state flip noise by restricting ourselves 
to only a subset of flips. One example would be a noisy 
channel with only 6 possible flips (k = 6) between two of 
the four basis states 
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and another a channel where only cyclic flips (k = 3) of 
all the basis states can appear 

A = /®i ; / 2 = / 2 ®i ; / 3 = / 3 ®i (31) 
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FIG. 10: (Color online) Quantum correlations between the x 
and y spatial dimensions for the Grover walk in the presence of 
a state-flip noise channel with p — 0.1. The noise is modelled 
as state-flips including all possible flips (k = 23), flips between 
only two of the basis states (k = 6) and cyclic flips of all the 
four basis state (k = 3) 





(a) 



(b) 



FIG. 9: (Color online) Quantum correlations between the 
x and y spatial dimensions created by the two-state walks 
for different bit-flip noise levels applied after evolution of one 
complete step. 



FIG. 11: (Color online) Probability distribution of (a) the 
Grover walk and (b) the two-state walks when subjected to 
a depolarizing channel. For the Grover walk the channel is 
given by Eq. (33) and for all walks the noise level is p = 0.1. 



The distribution is shown after 15 steps of evolution. 



with / 



The probability distributions for 



these two approximations are visually very similar to the 
situation where all possible state-flips are taken into ac- 
count (k = 23) and in Fig. 10 we compare the results 
obtained for the the x — y spatial quantum correlations 
for a noise level of p = 0.1. 

One can see that the spatial quantum correlations are 
affected stronger by the k = 3 than by the full k = 23 
flip noise. This implies that the two- and three state 
flips included in k = 23 acts as reversals of cyclic flips, 
thereby reducing the effect of noise. Since the trends 
for the decrease of the quantum correlation are function- 
ally similar for k = 3,6 and 23, we will use the model 
with k = 3 cyclic flips as the state-flip noise channel 
for the Grover walk in this section. Similarly, taking 
into account the computational limitations and following 
the model adopted for the state-flip, we will use a cyclic 
phase-flips to model the phase-flip noise, 



fi = r (g) t ; r2 = t 2 (g) i ; rs = f 3 (g) i, 



(32) 



where f 



with uj = e 4 



"1 " 

uj 
uj 2 
_0 u 3 _ 

Both these approximations for state-flip and phase-flip 
noise will makes the complex depolarization noise man- 
ageable for numerically treatment. The density matrix 
of the Grover walk can then be written as 



P_ 

15 

P_ 
15 



2=1 

3 



P_ 
15 



3 3 



i=l 3=1 

(l-p)[S 4 p 4s (t-l)Sj] 



1 3 



(33) 



and the pr obability distribution for this walk is shown in 
Fig. 11(a) for p = 0.1. The quantum correlations Q(p pp ) 
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FIG. 12: (Color online) Quantum correlations created by the 
Grover walk for different depolarizing noise levels. 



FIG. 13: (Color online) Particle-position quantum correla- 
tions created by the two-state walks for different depolarizing 
noise levels. 
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With increasing noise 



and Q(pxy) are shown in Fig. ' 
level, a decrease in Q(p pp ) is seen, whereas for the quan- 
tum correlations between the x and y spatial dimensions, 
Q(pxy), the same amount of noise mainly leads to a de- 



crease in the positive slope (see Fig. 12). From this we 
can conclude that the general trend in the quantum cor- 
relation due to state-flip noise (Fig.|4| and depolarizing 
noise is the same but the effect is slightly stronger when 
including the depolarizing channel. 



where &\ 

and Q(t - 1) 
d28bl. 
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and Q a (t — 1) are given by Eqs. (28a) and 



Similarly to the situation where we considered only 
state-flip noise after one complete step (see Figs. [8] and 
^) we find again that the quantum correlations for both 
walks behave nearly identical and only differ slightly in 
strength compared to the case of state-flip noise alone 
(see Figs. 13 and 14). 



2. Two-state walks 

The depolarizing channels for the alternate walk and 
Pauli walk can be written as 



hs(t) 



P2sa(t) 



I>Q(t-l)*J 
.<=i 

3 



+ (l-p)Q(t-l) 

(34a) 

+ {l-p)Qa{t-l), 
(34b) 



C. Robustness of two-state walk 

From the preceding sections we note that the x — y 
spatial correlations, Q(p xy ), have a larger absolute value 
for the two-state walks compared to the Grover walk and 
that the presence of noise affects all schemes in a similar 
manner. To quantify and better illustrate the effect the 
noise has we therefore calculate 



R(Pxy) 



Q(pxy) for noisy walk 



(35) 



Q(pxy) for noiseless walk' 
as a function of number of steps, which gives the rate 



of decrease in the quantum correlations. In Figs. 15 and 
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FIG. 14: (Color online) Quantum correlations between the x 
and y spatial dimensions created by the two-state walks for 
different depolarizing noise levels. 



16| we show this quantity in the presence of state-flip 
or depolarizing noise, respectively, for a noise level of 
p = 0.2. One can clearly see that in both cases the two- 
state walks are more robust against the noise at any point 
during the evolution. Note that the Grover walk only 
produces correlations from step 2 on, which is the reason 
for its graph starting later. 




FIG. 15: (Color online) Relative decay of the quantum cor- 
relations Q(pxy) in the presence of a state- flip noise channel 
with p 0.2. 




FIG. 16: (Color online) Relative decay of the quantum corre- 
lations Q(pxy) in the presence of a depolarizing noise channel 
with p — 0.2. 



V. BREAKDOWN OF STATE-FLIP AND 
PHASE-FLIP SYMMETRIES FOR FOUR-STATE 
WALKS 



The quantum walk of a two-state particle in ID is 
known to remain unaltered in the presence of unitary 
operations which equally effect each step of the evolu- 
tion. This is due to the existence of symmetries [52j [53] , 
which can help to identify different variants of the same 
quantum walk protocol and which can be useful in de- 
signing experimental implementation. For example, in 
a recent scheme used to implement a one-dimensional 
quantum walk using atoms in an optical lattice [56] , the 
conditional shift operator also flipped the state of the 
atom with every shift in position space. However, the 
existence of a bit-flip symmetry in the system allowed to 
implement the walk without the need for compensation 
of these bit-flips. In this section we look at possible sym- 
metries in the walks discussed above and show that the 
bit-flip and phase-flip symmetries, which are present in 
the evolution of the two-state particle are absent in the 
evolution of the four-state particle. 

The density matrix for a two-state quantum walk in 
the presence of a noisy channel will evolve through a 
linear combination of noisy operations on the state and 
the unaffected state itself. As an example we illustrate 
the symmetry due to bit-flip operations in the alternate 
walk with bit-flip noise after evolution of one complete 
step. The density matrix in this case is given by 



hs(t) 



&iS y S x p2s(t - l)SlSlo{ 

(1-P) \SyS X P 2s (t-l)StSl 



(36) 



where d\ 



1 

1 



) 1. When the noise level is p = 1 this 
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expression reduces to 



P2s(t) 



CTlSyS x p2s(t - 1)^5^5-} 

S%p 2s {t-l)Sl{S' y y. 



(37) 



where in the second line the bit-flip operation has been 
absorbed into the evolution operator S' y . This replaces 

|0)(0| and |1)(1| in S y by |1)(0| and |0Ul|, respectively. 
Similarly, for a phase-flip, (j\ in Eq. (37) is replaced by 

(T 3 = q ^ leading to |1)(1| in S y being replaced 

by — 11)(1| to construct S' y . 

An alternative way to look at this is to absorb the bit- 
flip or phase-flip operation into the coin operation. For a 
two-state walk using the Hadamard coin operation, the 
bit-flip operation after each steps corresponds to the coin 
operation taking the form, 



H' 



1 



y/2 I 



(38) 



and the phase-flip after each steps corresponds to the 
coin operation taking the form, 



H" 



1 

7i 



i i 
-i i 



(39) 



For a bit-flip or phase-flip of noise level p = 1 the 
Hadamard coin operation H can therefore be recast into 
a noiseless (p = 0) quantum walk evolution using H' 
and H" as coin operation. A noise level of p = 1 then 
returns a probability distribution equivalent to the noise- 
less evolution and consequently the maximum bit-flip and 
phase-flip noise level for a two state walk corresponds to 
p = 0.5. This is a symmetry within the alternate walk, 
which also holds for the Pauli walk. 

A state- flip noise channel for the four-state walk, on the 
other hand, evolves the state into a linear combination 
of all possible flips between the four basis states and an 
unchanged state for all values of p except for p = (see 
Eq.[24|). That is, only when p = 0, Eq. reduces to 



P4 s (t) = S 4 p4s(t - l)5 4 t , 



(40) 



whereas for any non-zero p including p = 1, Eq. (24) takes 
the form 



.1=1 



(41) 



Any attempt to absorb the noise operations fi into the 
shift operator or the coin operation leads to k differ- 
ent results, which have to be applied with probability |. 



Therefore, in contrast to the two-state evolution, a state 
flip noise level of p = 1 does not result in a pure evolu- 
tion equivalent to the situation for p = 0. However, if 
the state-flip noise is restricted to one possible operation 
(fi), the density matrix is no longer a linear combination 
of noisy operations and the unchanged state. Thus, for 
p = 1 in Eq. ( 24 ) a single noise operation can be absorbed 
into the GDO by changing the form of G (see Eq. (J5|). 

This absence of a useful symmetry for the four-state 
quantum walk reduces the chances to find an equivalent 
class of four-state quantum walk evolutions. Further- 
more, since the four-state quantum walk requires a spe- 
cific form of coin operation to implement the walk, any 
possible absorption will not result in a quantum walk 
in 2D, which is a significant difference to the two- state 
walks. 



VI. CONCLUSION 

In this work we have studied the decoherence prop- 
erties of three different schemes that realize a quantum 
walk in two-dimensions, namely the Grover (four-state) 
walk, the alternate walk and the Pauli walk. The noise 
for two-state particle evolution was modeled using a bit- 
flip channel and depolarizing channel. For the four-state 
evolution, different possible state-flip channels were ex- 
plored and we have shown a channel with 3 cyclic flips 
between all the four states can be used as a very good ap- 
proximation to the full situation. Similarly, we presented 
a possible model for the depolarizing channel of the four- 
state quantum walk. Using MID as a measure for the 
quantum correlations within the state, our studies have 
shown that the two-state quantum walk evolution is in 
general more robust against decoherence from state-flip 
and depolarizing noise channels. 

Following earlier studies on bit- and phase-flip symme- 
tries in two-state quantum walks in ID, we have shown 
that they also hold for two- state quantum walks in 2D, 
but break down for four-state 2D quantum walks. 

With the larger robustness against decoherence, the 
existence of symmetries which allow freedom of choice 
with respect to the initial state and the coin operation 
and the much easier experimental control, we conclude 
that two-state particles can be conveniently used to im- 
plement quantum walks in 2D compared to schemes using 
higher dimensional coins. An other important point to 
be noted is the straightforward extendability of the both 
two-state schemes to higher dimensions by successively 
carrying out the evolution in each dimension. 
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